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Abstract. In this article, we study finite dynamical systems defined over graphs, where the 
functions are applied asynchronously. Our goal is to quantify and understand stability of the dy- 
namics with respect to the update sequence, and to relate this to structural properties of the graph. 
We introduce and analyze three different notions of update sequence stability, each capturing dif- 
ferent aspects of the dynamics. When compared to each other, these stability concepts yield vastly 
different conclusions regarding the relationship between stability and graph structure, painting a 
more complete picture of update sequence stability. 



1. Introduction 

A graph dynamical system consists of (i) an underlying graph Y where each vertex has a 
state from a finite set K, [ii) a sequence of vertex functions that updates each vertex state as 
a function of its neighbors, and [iii) an update sequence that prescribes how to compose the 
functions to get the global dynamical system map. Examples of graph dynamical systems in the 



literature include systems with synchronous update such as cellular automata [|l8ll and Boolean 



networks y/ZO, those with asynchronous updates, e.g., sequential dynamical systems 1160 and 
functional linkage networks [9], and many probabilistic versions of the above systems. In this 
article, we focus on systems with asynchronous schedules, with the broad goal of understanding 
the sensitivity of the global dynamics with respect to changes in the update sequence. 

First, it is necessary to define what update order stability means, and how to quantify such a 
concept. There are many possible answers to this question, and there are vastly different ways to 
do this based on existing literature. After a brief introduction of the main concepts in Section [2l 
we devote the following three sections to presentations of three different aspects of update se- 
quence stability, each focusing on different properties of the dynamics. In Section [3l we look at 
the structure of the phase space and its robustness to changes in the update sequence. In Sec- 
tion m we analyze the number of attractor cycles that can be reached from a given initial state 
under different asynchronous update sequences. Finally, in Section [5l we look at the collection 
of periodic states and how they change (setwise) under different update sequences. 

Many of the concepts described above are scattered throughout the literature, usually studied 
outside of the theme of update sequence stability. This paper is the first of its kind to bring these 
ideas together with this common goal in mind. It is unrealistic to expect to develop a complete 
unifying theory of update order stability due to the sheer size, diversity, and inherent randomness 
of complex systems. However, there is still much be gained from such a study. One of the 
themes of this paper is understanding how certain stability properties are generally correlated 
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with physical properties of the underlying graph, such as edge density. For example, are systems 
over sparse graphs or dense graphs more stable? Not surprisingly, different notions of stability 
can lead to vastly different conclusions. We hope that this paper will clarify many of these ideas 
and set the groundwork for future investigations on the stability of asynchronous finite dynamical 
systems. 



2. Preliminaries 

Let Y be an undirected graph with vertex set v[F] = {1, . . . ,n}. Each vertex i has a state 
Xi E K, where i^' is a finite set (frequently K = ¥2 = {0, 1}), and a vertex function fi that 
updates its state as a function of the states of its neighbors (itself included). When applying these 
functions asynchronously, it is convenient to encode fi as a Y -local function Fi : K'"- — K^, 
where 

In the remainder of this paper, we will focus on a class of graph dynamical systems with asyn- 



chronous update schedules, called sequential dynamical systems alw . These systems have a 
prescribed update schedule dictated by a sequence of vertices cu G Wy, where Wy is the set of 
non-empty words over v[Y]. 

Definition 2.1. A sequential dynamical system (SDS) is a triple {Y, ^y,uj), where Y is an undi- 
rected graph, ^y a sequence of F-local functions, and uj e Wy a length-m update sequence. 
The SDS map is the composition 

[dy, Uj] = o F^^_^ o ■ ■ . o o . 

Henceforth, we will primarily restrict our attention to the case where to E Sy, the set of 
permutations of v[Y], i.e., words where each vertex appears precisely once. To emphasize a 
permutation update sequence, we will use vr instead of to. All of the concepts carry over to 
general asynchronous systems under slight modification, though not much additional insight is 
gained in doing so, whereas the notation becomes more involved and the main concepts obscured. 

The phase space of a finite dynamical system 0: — > K'^ is the directed graph r(0) with 
vertex set and edges (x, for each x E K". Clearly, r(0) consists of disjoint cycles 
(periodic states), with directed trees (transient states) attached. 



3. Equivalence of Dynamics 

Any reasonable notion of stability of dynamics must be coupled with a notion of equivalence 
of dynamics. For dynamical systems, a natural place to begin is with the structure of the phase 
space, since this encodes the global dynamics. We present here several different types of equiv- 
alence, which vary in what is meant by having equivalent structure. 

Definition 3.1. Let 0, ip : if" — > K'^ be finite dynamical systems. Then and ^ are: 

(1) functionally equivalent if (p = ip. 

(2) dynamically equivalent if (f) o h = h o ip for some bijection h : — > 

(3) cycle equivalent if 0|per(0) o h = ho '0|per(^) for some bijection h: Per(0) — > Per('0). 
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Clearly, cf) and ip are functionally equivalent iff their phase spaces are identical, dynamically 
equivalent iff their phase spaces are isomorphic, and cycle equivalent iff their phase spaces are 
isomorphic when restricted to the periodic states. For a fixed sequence '^y of F-local functions, 
different update sequences give rise to different SDS maps, and each notion of equivalence men- 
tioned above partitions the set of all such maps into equivalence classes. The number of SDS 
maps up to equivalence obtainable by varying the update sequence is a measure of stability. As 
we will see, each type of equivalence of dynamics corresponds with a combinatorial equivalence 
relation on the set of update sequences Sy (or Wy), which in turn is strongly tied to the struc- 
ture of the underlying graph Y. Thus, knowledge of the structure of these equivalence classes is 
central to quantifying update sequence stability, and its relationship to the base graph. We will 
summarize this below. 

3.1. Functional Equivalence. Every update sequence n e Sy defines a partial ordering <^ of 
v[Y], and this is naturally represented via the set Acyc(F) of acyclic orientations on Y. Explicitly, 
TT gives rise to the acyclic orientation 0-,^ of Y obtained by orienting each edge {i,j} of Y as 
(i, j) if i occurs before j in tt and as (j, i) otherwise. The partial order <t, is defined by i <t, j iff 
there is a directed path from i to j in Oy. If vr, a G Sy are linear extensions of the same partial 
ordering, then [^y, vr] = ['^y, a], and we say that vr ~q, a. In fact, we have a bijection 

(pa-. Sy/r^^ — . Acyc(r) , 

hence a{Y) := | Acyc(F) | is an upper bound for the number of distinct (permutation) SDS maps 
over ^y up to functional equivalence. This bound is sharp (see and it is counted by the 
Tutte polynomial, as ol{Y) = Ty(2, 0). Additionally, the bijection (pa allows us to construct a 
transversal of update sequences from Sy, rather than considering all n\ update sequences. 

3.2. Dynamical Equivalence. Let Aut(F) denote the automorphism group of Y, and define the 
^n-action on K"- by a{xi, . . . , Xn) = {x^-i(^i), . . . , x^-ii^n))- The group Aut(F) acts on the set 
Sy/^a by 7 ■ [Tr]a = [yTc]a- If dv is Aut(F)-invariant (i.e., 7 o F^, = F^i^y) o 7 for each local 
function, and each 7 G Aut(F)), then two SDS maps over update sequences in the same orbit are 



conjugate, by 70 [^y, vr] 07 ^ = [^y , 77r] (see nlw ). This puts an equivalence relation on Sy (and 



hence on Acyc(F)) that we denote by Since conjugate maps are dynamically equivalent, the 
function aiY) := |Acyc(F)/ ~a | is an upper bound for the number of SDS maps over ^y up 
to dynamical equivalence. This quantity, which counts the number of orbits in Acyc(F) under 
Aut(F), can be computed though Burnside's lemma, as 

Here, (7) \ Y denotes the orbit graph of the cyclic group G= (7)andF, (see 10, i). We believe 
this bound is sharp, and have shown this for special graph classes, but have no general proof. 

3.3. Cycle Equivalence. To see how cycle equivalence arises in the SDS setting, we first define 
a way to modify an acyclic orientation, called a source-to-sink operation, or a click. This consists 
of choosing a source vertex and reversing the orientations of all incident edges. The reflexive 
transitive closure of this puts an equivalence relation on Acyc(F) (and on Sy), denoted ~k- Thus, 
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two acyclic orientations related by a sequence of clicks are K-equivalent. In [|15l1 . it is shown that 
if TT ~K a, then [^'y, tt] and [5'y, cr] are cycle equivalent. Therefore, k{Y) := |Acyc(F)/~K I is 
an upper bound for the number of SDS maps over up to cycle equivalence. This quantity is 
also counted by the Tutte polynomial, by k{Y) = Ty(l, 0) (see [14]). As with a-equivalence, 
there is a nice transversal for /t-equivalence, arising from the bijection 

(3.1) 0.,.: Acyc„(r) ^ Acyc(F)/~« . 

where Acyc^(F) is the set of acyclic orientations of Y where vertex v E v[Y] is the unique 
source. 

If additionally, Aut(F) is nontrivial, we can use the action of this group on Acyc(y) / ~k to 
impose a coarser equivalence relation ~s on Acyc(F). This is clear, since dynamically equivalent 
SDS maps are trivially cycle equivalent. The quantity k{Y) := |Acyc(F)/ ~k |, which is the 
number of orbits in Acyc(F)/ ~k under Aut(F), is therefore a stronger upper bound for the 
number of SDS maps over up to cycle equivalence when ^y is Aut(y)-invariant. 

3.4. Summary of Equivalence. To summarize, for each of the following notions of equivalence 
of SDS maps, there is a corresponding equivalence relation on Sy (and hence on Acyc(F)): 

■ Functional equivalence. Relation: ~a with transversal Acyc(F) via^Q,; 

■ Dynamical equivalence. Relation: with transversal Acyc(F)/~a; 

■ Cycle equivalence. Relation: ~k and transversal Acyc^{Y) via (p^^^ for v E v[Y]. 

In each of these three cases, update sequences that are equivalent give rise to SDS maps with 
equivalent dynamics. Additionally, it is clear that functionally equivalent systems are dynami- 
cally equivalent, which in turn are cycle equivalent. In the general case where the functions ^y 
are not Aut(y)-invariant, or when Aut(F) = 1 (which is true of almost all random graphs [l6|]), 
a-equivalence reduces down to a-equivalence. The difference between a- and a-equivalence 
only captures special symmetries in the graph. In light of this, we will focus on a(Y) and k{Y) 
for our discussion of stability, though «(F) was mentioned for sake of completeness. 

Since the functions a(Y) and k(Y) bound the number of SDS maps obtainable up to equiva- 
lence under different update sequences, they can be considered as measures of stability. Comput- 
ing these quantities is in general, NP-hard, since they are evaluations of the Tutte polynomial {t\\. 
However, even though we can not compute these functions exactly for large networks, an un- 
derstanding of the mathematics involved provides valuable insight about update order stability. 
For example, consider two extreme cases: (i) Y is an n-vertex tree, and (ii) Y = Kn, the 
complete graph on n vertices. In the first case, k{Y) = 1, and in the latter, k(Y) = (n — 
Moreover, since a and k are Tutte- Grothendieck invariants, if Z < Y then a(Z) < a(Y) and 
k(Z) < k{Y), i.e., a and k monotonically increase as additional edges are introduced to the base 
graph. Two intermediate cases (between the tree and complete graph) were studied in iflsll - the 
graphs corresponding to the radius-1 and radius-2 cellular automata rules. For radius-1 rules 
(functions over the the circular graph Circ„), there are at most k{Y) = 0{n) cycle structures for 
(permutation) SDS maps. In contrast, an SDS with a radius-2 rule has base graph Y = CirCn 2, 
the circular graph where every vertex is additionally connected to its distance-2 neighbors. In 
this case, k(Y) = 0{n ■ 2"). Though we cannot compute k.{Y) explicitly for most graphs, we 
see a clear correlation between edge density of the base graph and the stability of a dynamical 
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system over it. Thus, using a or ^-equivalence as a notion of stability, we can make the general 
statement: 

Update order stability and edge density are positively correlated. 

3.5. An Example. We conclude this section with an example illustrating the difference between 
the various notions of equivalence described in this section. 

Example 3.2. Let Y = Circ4, with vertex functions nors: F2 — > ¥2 where nor3(0,0,0) = 1 
and nors{x,y, z) = 1 for all other arguments. Using update sequences tt = (1,2,3,4), n' = 
(1, 4, 2, 3) and n" = (1, 3, 2, 4), the phase spaces are shown in Figure [T] None of these SDS 




Figure 1. The phase spaces of Example 13. 2[ 

maps are functionally or dynamically equivalent, but [Nory, tt'] and [Nory, tt"] are cycle equiv- 
alent, as illustrated in the two rightmost phase spaces in Figure [T] In fact, it is elementary to 
compute K(Circ4) = 2, thus the two cycle configurations in Figure [T] are the only possible such 
configurations up to isomorphism when Y = Circ4 and K = ¥2. 

4. Limit Sets and Reachability 

Another natural measure for update sequence stability is a quantification of how many attractor 
cycles can be reached from a certain initial state. There are many ways to develop such a measure 
- one can compute the average number of periodic states reachable from a random starting point, 
the maximum number of reachable periodic states, the number of reachable attractor cycles, pos- 
sibly normalized for size, etc. For general complex networks, such computations are intractable, 
and Monte Carlo type methods may be an alternative. However, it is at times possible to esti- 
mate or bound such quantities asymptotically. In this section, we will summarize some recent 
results of one these measures that was developed to study a the robustness of a gene annotation 
algorithm that basically uses a threshold SDS [9]. 

Threshold functions are common in discrete dynamical system models fid] . A Boolean thresh- 
old function on k arguments is a function fk,rri - ¥2 — > ¥2 such that f{xi, . . . ,Xk) = 1 when 
at least m of the arguments Xi are 1, and f{xi, . . . , Xk) = otherwise. An SDS is a thresh- 
old SDS if every vertex function is a threshold function. It is well-known that limit sets of 
threshold SDSs contain only fixed points llldl . and additionally, that fixed points are independent 
of update sequence. This class plays a special role in complexity theory since in many cases. 
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threshold functions are the most general function class for which problems are computationally 
tractable 

Functional linkage networks (FLNs) and the gene annotation algorithm of ^ were among 
the motivations for considering the type of update sequence stability measures discussed in this 
section. An FLN (for a fixed gene ontology function /) is a graph where the vertices represent 
proteins, and an edge is present if two proteins are thought to share the same function, with 
an edge weight Wij describing certainty. The vertices are assigned states from K = {1,-1,0} 
denoting whether they are annotated with / (state 1), are not annotated with / (state —1), or if 
annotation is unknown (state 0). Each vertex i is updated by changing its state to 1 if ^ WijXj > 
T, where the sum is taken over all neighboring vertices, and T is some fixed threshold, and to 
— 1 otherwise. The algorithm updates the vertices in some fixed order chosen at random, and 
terminates when a fixed point is reached. 

The FLN algorithm uses an asynchronous update schedule to ensure that a fixed point is 
reached (there may be periodic cycles of length > 1 under a synchronous update schedule). 
In light of this convenient choice of update scheme, and the random choice of update sequence, 
one would hope that such a gene annotation algorithm would exhibit update order stability. In 
other words, the final (fixed point) state reached should be robust with respect to changes in up- 
date sequence and perturbations of the initial state, in that the fixed point reach has little or no 
dependence on the update sequence. Unfortunately, such stability cannot always be guaranteed. 
As an example, there are general classes of threshold SDSs sharing many of the essential FLN 
algorithm properties, but which exhibit exponential update sequence instability [11] for many 
choices of the initial state. Not surprisingly, such instability strongly depends on the choice of 
underlying graph. 

The co-limit set of 0: — > from x E is the set to {(f); x) of periodic states x E K"^ 
such that (f)'^(x) = y for some m > 0. This is a specialization of the classical definition of 
cu-limit set to the case of finite phase spaces. If V C Sy is a collection of update sequences, then 
the a;-limit set of from x with respect to V is defined as 



If there are large periodic cycles in the phase space of the SDS, then the size of an cj-limit set 
may not necessarily be too insightful. However, if there only fixed points, as is common in many 
applications, then u;p(s) is a direct measure of stability that describes how many fixed points can 
be reached from x by choosing the update schedule from V. For a set of functions ^y, define 



This counts the maximum possible number of periodic points that can be reached from any state 
by variation of the update order. As in the previous section, we begin by considering two extreme 
cases; when the base graph is a tree, and when it is the complete graph Kn. The following result 



Theorem 4.1. Let ^y be a sequence of 2-threshold functions. If Y is the star-graph on n vertices, 
then uidv) = 2" - n. If F = K^, then loCSy) =n + l. 




uoidv) = max\{ujsy{x) \ x E fsT"}] . 



is proven in [iTlll . 
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Additionally, this is extended to the random graph model G{n,p) for ranges of p relevant to 
many applications. The main result is that for sparse graphs, uj{dY) = 0(2") (with high proba- 
bility), but for dense graphs, uj{^y) = Q{n). While these systems are deliberately constructed 
to possess dynamical instability, they nonetheless offer valuable insight into possible dynamics 
of such systems. With this notion of stability with respect to update sequence, one can make the 
following general statement, which is in direct contrast to the conclusion of the previous section: 



Update order stability and edge density are negatively correlated. 



5. PERIODIC Points and Word Independence 

We come at last to our third and final stability concept for asynchronous dynamical systems, 
which stems from recent work on word independent systems [js, 12, 13]. Like cycle equivalence. 



this considers invariance properties of the periodic points of SDSs, rather than the transient states. 
A sequence of F-local functions is word independent if the set of periodic states Per([5^y , uj]) 
is the same for all complete words u G Wy (words over v[F] where each vertex occurs at least 
once). Note that word independence only considers the periodic states as a set, and ignores the 
orbit structure. Word independence was first studied in [8|], where it was shown that most clas- 
sic symmetric functions such as the logical AND, OR, NAND, NOR functions, along with their 
sums, were word independent, for asynchronous cellular automata (i.e., over the circular graph). 
Other common functions, such as parity, minority, and majority functions had this property as 
well. Recently, in a series of two papers 13], we have expanded this study to all 256 elemen- 



tary cellular automata rules, and classified the word independent ones, which are independent of 
the size of the underlying (circular) graph. 

Theorem 5.1. Precisely 104 of the 256 elementary cellular automata rules induce a word inde- 
pendent SDS over Circ„, for all n > 3. 



The classification of these 104 rules is contained in [1131] . and the structure is further developed 
in lfl2|] with an analysis of their dynamics groups, which describe how the periodic states are 
permuted within these systems as the local functions are applied. The latter paper additionally 
contains more general results on dynamics groups. For example, when the state space h K = 
{0, 1}, as is typical in finite dynamical systems research, the dynamics group is the homomorphic 
image of a Coxeter group and for general K, an Artin group. 

The papers above are only the beginning of this line of work. There is much to be done to 
better understand systems that fail to be word independent. Invariance of periodic states under 
different update sequences can be used as a notion of update sequence stability. This can be 
further quantified as follows. When a sequence of functions is not word independent, we can 
define a measure that captures how "close" it is to being word independent, by defining 



X e Per[g'y, Tc] for all n e Sy\ 



I X G Per[^Y, tt] for some tt E Sy] 

Clearly, < pidr) < 1, with p(^^y) = 1 iff is word independent. This quantity measures 
how stable the periodic point sets are as a whole, to changes in the update sequence. Under 



8 



MATTHEW MACAULEY AND HENNING S. MORTVEIT 



this definition, word independent sequences of functions are the most stable. In general, it is 
difficult to determine if a sequence of functions is word independent. Though there may exist a 
correlation between p(5^y) and the topology (e.g., edge density) of Y, there is no good reason to 
believe that such a link should exist. Rather, there are many examples of functions that are word 
independent for all graphs, such as threshold functions, monotone functions, functions with only 
fixed points, parity functions, and the logical NOR functions. Therefore, it is not far-fetched to 
make the following hypothesis: 

There is little, if any, correlation between update order stability and edge density. 

6. Concluding Remarks 

This article should raise a red flag for anyone anyone attempting to study and characterize 
stability with respect to update sequence variation. We saw how different notions of stability can 
lead to different conclusions of the relationship to graph structure. In summary, we looked at 
three properties of the dynamics and how robust they are to changes in the update sequence. In 
particular, we considered: 

■ The variation of structural properties of the phase space as the update sequence changes; 

■ The variation of reachable limit sets as the update sequence changes; 

■ The variation of the sets of periodic states as the update sequence changes. 

For each of these notions, we defined natural measures quantifying stability, and then inquired 
about the relationship between these measures and the structure of the underlying graph. Though 
such a question is quite open-ended, we began by looking at the extreme cases of sparse graphs 
(e.g., trees) and dense graphs (e.g., complete graphs) to gain some insight, before considering 
intermediate cases. Up until now, these three notions of stability have only been studied inde- 
pendently, and moreover, outside of the setting of stability. In this article, we see first-hand how 
differently they behave with respect to a basic property such as edge density of the base graph. 
Perhaps there are ways to draw connections between them to paint a more complete picture of 
stability. It is our hope that this article will help lay the groundwork for such future research. 
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